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Abstract 

Polydnaviruses are obligate symbionts integrated as proviruses in the genome of some ichneumonoid wasps that parasitize 
lepidopteran larvae. Polydnavirus free viral particles, which are injected into the host at oviposition, express virulence factors 
that impair immunity and development. To date, most studies have focused on the molecular mechanisms underpinning 
immunosuppression, whereas how viral genes disrupt the endocrine balance remains largely uninvestigated. Using 
Drosophila as a model system, the present report analyzes the function of a member of the ankyrin gene family of the 
bracovirus associated with Toxoneuron nigriceps, a larval parasitoid of the noctuid moth Heliothis virescens. We found that 
the TnBVankl expression in the Drosophila prothoracic gland blocks the larval-pupal molt. This phenotype can be rescued 
by feeding the larvae with 20-hydroxyecdysone. The localization of the TnBVANKl is restricted to the cytoplasm where it 
interacts with Hrs and Alix marked endosomes. Collectively, our data demonstrate that the TdBVANKI protein acts as a 
virulence factor that causes the disruption of ecdysone biosynthesis and developmental arrest by impairing the vesicular 
traffic of ecdysteroid precursors in the prothoracic gland steroidogenic cells. 
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introduction 

Parasitic wasps represent tiie largest group of parasitoid insects 
which attack and parasitize a number of insect species, exploiting 
different developmental stages [1]. These parasitic insects have a 
peculiar injection device, the ovipositor, which is used to deliver 
the egg along with host regulation factors that primarily disrupt 
the host immune reaction and endocrine balance to create a 
suitable environment for the development of their progeny [2,3]. 
These host regulation factors include viruses of the Polydnaviridae 
family, obligate symbionts of ichneumonid and braconid wasps 
attacking larval stages of lepidopteran hosts, and respectively 
classified in the genera Ichnovims (TV) and Bracovirus (BV) [4] . 

Polydnaviruses (PDVs) [5,6] are integrated as proviruses in the 
genome of parasitoid wasps and their transmission to offspring is 
strictiy vertical, through the germline. The genome encapsidated 
in the viral particles is made of multiple circular dsDNA segments, 
which have an aggregate size ranging between 190 and 600 kb. 
PDVs only replicate in the epithelial cells of the calyx, a specific 
region of the ovary, where they accumulate to a high density to be 
injected at oviposition along with the venom and the egg. 
Free PDV particles infect the host tissues without undergoing 



replication, and express virulence factors that alter host physiology 
in ways essential for offspring survival [5]. 

Evolutionary convergence of independent host-virus associa- 
tions has favored the selection of gene families shared by both IV 
and BV [7-9]. For example, protein tyrosine phosphatases (PTP) 
and ankyrin motif proteins (ANK) are widely distributed in many 
PDVs and expressed to different degrees in virtually all host tissues 
analyzed so far, indicating that they play a key role in successful 
parasitism [10]. 

We know more about the functional bases underlying 
the immune disguise in parasitized hosts than we do about 
how the host developmental alteration is induced [3,1 1]. This 
is due to the complexity of the developmental mechanisms and 
to the concurrent action of various virulence factors, which 
often have redundant and overlapping effects on the regulating 
gene networks [2,3,12]. 

One of the best characterized developmental syndromes has 
been described in the host-parasitoid association Heliothis virescens- 
Toxoneuron nigriceps (Lepidoptera, Noctuidae - Hymenoptera, 
Braconidae) [13]. Briefly, in this experimental model the last 
instar larvae fail to pupate and show a higher nutritional suitability 
for parasitoid larvae. The developmental arrest is partly due to a 
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marked depression of ecdysone (E) biosynthesis by the prothoracic 
gland (PG), induced by the infection of the bracovirus associated 
with T. nigriceps (ThBV). The inhibition of E biosynthesis is further 
reinforced by the conversion of the very low amounts of 20- 
hydroxyecdysone (2()E) produced to inactive polar metabolites, a 
transformation mediated by teratocytes, special cells deriving from 
the parasitoid's embryonic membrane. 

The active transcription of TnSV genes in the PG of 
parasitized tobacco budworm larvae is required to disrupt their 
ecdysteroid biosynthesis, which remains very low and fails to 
increase in response to prothoracii'otropie hormone (PTTH) 
stimulation [13,14]. Unraveling the functional role of a specific 
virulence factor at molecular level is not easy when the natural 
host is used for these studies, due to the limited availability of 
genomic information and molecular tools. Therefore we used 
Dro.sophila melanogaster as an ideal experimental model to study 
TnBV genes. With this approach, we started the functional 
characterization of a member of the viral ankyrin [ank) gene family 
of TnBY, TnBVankI [15], showing that the expression of this gene 
in Drosophila germ cells alters the microtubule network function in 
the oocyte [16]. In the present study we analyze the effect of 
TnBVankI gene expression during Drosophila development. Inter- 
estingly, we found that TnBVankI expression in the PG cells 
blocks the transition from larval to pupal stage, mimicking the 
developmental arrest observed in H. virescens larvae parasitized by 
T. nigriceps. 

Materials and Methods 

Fly strains 

Stocks were raised on standard cornmeal/yeast/agar medium 
at 21°C and crosses were made at 25 °C unless otherwise stated, 
jjif/'^'^"' was used as the wild-type stock in this study. The UASp- 
TnBVankl strain (genotype: UASp-TnBVankl / UASp-TnBVan- 
kl ;UASp-TnBVankl / UASp-TnBVankl ; +/+) was generated in our 
laboratory [16]. 

The following stocks were obtained from the Bloomington Stock 
Center: hairy-Gal4 (#1734: w'; P{GawB}h'J''), UASp-a-tubulin- 
GFP (#7374: y' w*; P{UASp-GFPS65C-QtTub84B} 14-6-II), 
UAS-p35 (#5073: w*; P{UAS-p35.H}BH2) and tub-GalSO"" 
(#7019: w*; Plw"^""^ = tubP-GAL80''}20; TM2/TM6B, Tb). 

phantom-Gal4 and P0206-Gal4 were a gift from C. Mirth 
(phantom-Gal4, UAS-mCD8::GFP/ TM6B and P0206-Gal4, UAS- 
mCDSr.GFP), august21-Gal4;phantom-Gal4 was kindly provided by 
M. Jindra (aiigust21-Gal4/ CyO;phantx)m-Gal4/phantom-Gal4). 

The stocks used for Gal4 driven expression of UASp-TnBVankl 
referred in Figure S2 and listed in Table SI are from Bloomington 
Stock Center. 

Crosses 

Females UASp-TnBVankl were crossed to males of the different 
Gal4 lines. As control, females jw''^'^' were crossed to males of the 
same Gal4 lines. 

For microtubules analysis, females UASp-TnBVankl were 
crossed to males UASp-a-tubulin-GFP; phantom-Gal4. 

For GalSO*^ experiment, females UASp-TnBVankl were 
crossed to males tub-Gal80'';phantom-Gal4,UAS-mCD8::GFP/ 
TM6B and females yw*"''^'^ were crossed to males tub-Gal80'" 
';phantom-Gal4,UAS-mCD8::GFP/TM6B as control. 

To coexpress p35 and TwBVANKl in PG cells females UASp- 
TnBVankl ;UASp-TnBVarJcl;UAS-p35 were crossed to males phan- 
tom-Gal4, UAS-mCD8::GFP/ TM6B. 



Larval length measurements 

Five UASp-TnBVankl /+;UASp-TnBVankl /+;hairy-Gal4/+ larvae 
at different days after egg deposition (AED) and five control larvae 
were ice-anesthetized and photographed using a Nikon Eclipse 90i 
microscope. Images were taken at 4X magnification and the larval 
length was measured with NIS-Elements Advanced Research 3. 10 
software. 

20E titer 

Five larvae at different developmental stages were collected and 
washed with PBS and immediately frozen by liquid nitrogen. 
Samples were added 200 nl of methanol, homogenized and 
transferred into 1.5 ml plastic tubes. After 10 minutes centrifuga- 

tion (12,000 rpm at 4°C) the supernatant was collected into a new 
tube, the precipitate was re-extracted with 200 |J.l of methanol and 
the supernatant was added to the previous one. After 30 minutes 
on ice, the samples were centrifuged following the same 
conditions. Samples were dried to remove methanol and then 
dissolved in the borate buffer. The standard cur\'e was generated 
according to the standard process of the RIA protocol [1 7] and 
then the 20E titer in samples was calculated. 

Rescue experiment 

UASp- TnB Vankl /+; UASp- TnB Vankl /+;phantom-Gal4, UAS- 
mCD8.:GFP/+ and control larvae were collected at 106 h AED and 
placed in three groups of ten individuals at 25°C in new tubes 
supplemented with 20E (Sigma) dissolved in ethanol at 1 mg/ml. 
Control larvae were fed only with ethanol. 

Prothoracic gland and cellular size measurements 

For measurements of PG area and its cellular size, confocal 
images of 50 PGs taken at 40X magnification were quantified with 
ImageJ software. 

Statistical analysis 

Statistical comparison of mean values was performed by 
unpaired t-test, using GraphPad Prism 4 software. 

Immunofluorescence microscopy 

Larvae were dissected at room temperature in IxPBS pH 7.5 
(PBS) and fixed in 4% formaldehyde for 20 minutes at room 
temperature. After three washes in PBS, larvae were permeabi- 
lized in PBT (PBS pH 7.5-1-0.3% Triton X-100) for 1 h, washed 
three times 5 minutes each in PBT and 10 minutes in PBT-H 
2%BSA solution. After that, the larvae were incubated, overnight 
at 4°C, with primary antibodies diluted in PBT-h2%BSA. Larvae 
were washed three times 10 minutes each in PBT, 10 minutes in 
PBT-I-1%BSA solution and incubated 2 hours at room tempera- 
ture on a rotating wheel with secondary antibodies diluted in 
PBT-Hl%BSA. After several washes in PBT, the ring glands were 
dissected and mounted on microscopy slides in Fluoromount G 
(Electron Microscopy Sciences). Subsequendy samples were 
analyzed by conventional epifluorescence with a Nikon Eclipse 
901 microscope or with TCS SL Leica confocal system. Images 
were processed using Adobe Photoshop CS4 and Adobe Illustrator 
CS4. 

TRITC-Phalloidin staining was carried out, after incubation 
with secondary antibodies, by washing larvae three times with PBS 
and then by incubating larvae for 20 minutes with TRITC- 
PhaUoidin (40 Hg/ml in PBS, Sigma). 

For Propidium Iodide nuclear counterstaining, the larvae were 
treated with RNase A (400 Hg/ml in PBT, Sigma) overnight at 
4°C. After three washes in PBT, the larvae were labeled for 2 
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hours with Propidium Iodide (10 Hg/ml in PBT, Molecular 
Probes). 

The following primary antibodies were used: polyclonal rabbit 
anti-Dib 1:200 [18], anti-Cleaved Caspase-3 1:25 (9661, Cell 
signaling Technology), anti-Rab7 1:2000 [19] and anti-Rabll 
1:5000 [19] were detected with DyLight 649-conjugated goat anti- 
rabbit 1:500 Jackson). Polyclonal rabbit anti-TwBVANKl 1:200 
[16] was detected using Cy3- (1:1000) and DyLight 649- (1:500) 
conjugated goat anti-rabbit (Jackson). Monoclonal mouse P1H4 
anti-Dynein heavy chain 1:200 [20], anti-Rab5 1:25 (610281, BD 
Biosciences) and anti-Alix 1:100 [21] were detected with Cy3- 
conjugated goat anti-mouse 1:1000 (fackson). Polyclonal guinea 
pig anti-Hrs 1:1000 [22] was detected with DyLight 649- 
conjugated goat anti-guinea pig 1:500 (Jackson). 

Terminal deoxynucleotidyl transferase-mediated dUTP 
Nick End Labeling (TUNEL) 

Five days AED larvae were dissected at room temperature in 
PBS, fixed in 4% formaldehyde for 20 minutes according to a 
protocol previously described [23]. After TUNEL incubation, 
anti-Digoxigenin 1 : 1 00 (Roche) was detected with Cy3-conjugated 
goat anti-mouse 1:1000 (Jackson). 

Filipin and Oil Red O staining 

Ring glands were fixed in 4% formaldehyde for 20 minutes and 
washed three times in PBS for 5 minutes each. Samples were 
stained with 50 |xg/ml of filipin (Sigma) for 1 h or incubated in an 
Oil Red O (Sigma) solution at 0.06% for 30 minutes. After 
incubation tissues were washed twice with PBS before mounting in 
Fluoromount-G. Samples were analyzed by conventional epifluor- 
escence with a Nikon Eclipse 90i microscope or with a Nikon 
Eclipse 901 confocal microscope. Images were processed using 
Adobe Photoshop CS4 and Adobe Illustrator CS4. 

Colocalization analysis 

Thresholds of confocal images were set in Adobe Photoshop 
CS4 to exclude background staining. 509 Hrs positive vesicles 
were analyzed per TnBVANKl and Hrs staining. 443 
TnBVANKl positive vesicles were analyzed per TnBVANKl 
and Alix staining. 1 1 8 Hrs positive vesicles were analyzed per Alix 
and Hrs staining. 

Images were processed with the CDA plugin of ImageJ to 
obtain Pearson's coefficient (from -1-1 = complete correlation, to — 
1 = anti-correlation with 0 = no correlation) [24] . 

Results 

Expression of TnBVankI in the prothoracic gland induces 
developmental arrest at third instar larvae 

TnBVankI gene expression during Drosophila development was 
targeted with the GAL4/UAS binary system [25]. We used a 
transgenic Drosophik stock carrying two copies of the TnBVankI 
gene under the control of the UASp sequences [16]. Expression of 
this transgene was induced using different Gal4 drivers. Our first 
analysis expressed the TnBVankI transgene during (-mbrv onic and 
larval development using the hahy-GaM driver [h-Gal4) [25]. 
TnBVankI expression did not appear to affect embryonic and 
larval development but, interestingly, all larvae failed to pupate 
and died after an extended third instar larval life, which lasted up 
to three weeks (Figure 1 A). By measuring larval size, we found that 
four days AED the larvae expressing TnBVankI did not 
significantly differ from control yw;h-Gal4 (n = 5; t = 0.8557; NS) 
(Figure IB). Moreover, they continued to feed and significantly 



increased in size during their prolonged larval life, reaching at 
eighteen days the maximal length (Figure 1A,C; n = 5; t = 6.765; 
p<0.0001), while control regularly pupated on day six AED 
(Figure lA). Since h-Gal4 is expressed in various larval tissues, the 
observed developmental arrest suggested us that TnBVankI 
expression could ha\ e reasonably affected the ring gland function, 
the major site of production and release of developmental 
hormones. 

The Drosophila ring gland (Figure 2A) consists of the prothoracic 
gland (PG), which is composed of steroidogenic ( ells synthesizing 
the E, the corpora allata (CA) that produce- the juvenile hormone, 
and the corpora cardiaca (CC), which play a key role in the 
regulation of metabolic homeostasis [26]. We assessed if the 
targeted expression of the TnBVankI gene using different ring 
gland Gal4 drivers (Figure 2B) was able to reproduce the effect 
observed when the transgene was expressed using the h-Gal4. 
When the TnBVankI gene was expressed in both CA and PG cells, 
using the P0206-Gal4 driver, all the larvae failed to pupate and 
showed the same phenotype obtained with h-Gal4. Conversely, 
when the august2TGal4 {aug2TGal4) driver specifically targeted the 
expression of TnBVankI in the CA, no effects on developmental 
timing were observed and regular progeny were obtained. 
Moreover, we specifically induced expression of the TnBVankI 
gene in the PG using the phantom-Gal4 {phm-Gal4) driver, which is 
strongly expressed in this gland. None of the larvae pupated and 
they had an extended larval Ufe as shown using the P0206-Gal4 
driver (Figure 2B). These data indicate TnBVANKl impairs PG 
function causing the block of larval-pupal transition. 

We specifically expressed TnBVankI in several other tissues, 
using different Gal4 drivers, and monitored the timing of 
development and the adult phenotype, which were in all cases 
not affected (Table SI). 

Collectively, our data suggest that the expression of 
TnBVANKl has the potential to interfere with the steroid 
biosynthesis, as further indicated by the targeted expression of 
this viral ANK protein in PG, which is characterized by 
developmental arrest of mature larvae and the absence of a 
systemic injury response [27]. 

Ecdysteroid biosynthesis is impaired in TnBVankI larvae 

To assess whether the developmental arrest induced by 
TnBVANKl was due to a reduced level of 20E, we measured 
the whole body 20E titer in larvae expressing TnBVankI by the 
phm-Gal4 driver and in control larvae (Figure 2C). At 1 10 h AED, 
wild type third instar larvae enter the wandering stage and, at 
25°C, they become white pre-pupae at 120 h, after the surge of a 
20E peak [17]. At 120 h AED and during their abnormal 
extended larval life, the 20E levels measured in phm-Gal4> 
TnBVankI larvae are extremely reduced and significantly lower 
than that measured both in UASp-TnBVankl larvae (n = 5; 
t= 10.12; p<0.0001) and in phm-Gal4/TM6B larvae (n = 5; 
t = 8.196; p<0.0001). 

To further demonstrate that the block of the transition to pupal 
stage showed by the phm>TnBVankl larvae (hereinafter TnBVankI 
larvae) was actually due to a low level of 20E, we carried out a 
20E-feeding rescue experiment. Third instar TnBVankI larvae 
were fed with yeast paste containing 2()E dissoh-c'd in ethanol at 
106 h AED, just before the onset of the ecdysteroid peak occurring 
in the wild-type. As expected, at 120 h AED, 1Q% of control larvae 
started to pupate and within the following 20 h all of them reached 
the pupal stage (n = 30). Pupation of TnBVankI larvae fed with 20E 
followed almost an identical pattern, with 100% pupation (n = 30) 
attained only 1 day later, but failed to progress to the pharate stage 
(Figure 2D). Instead, TnBVankI larvae treated only with yeast and 
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h-Gal4>TnBVank1 






Figure 1 . TnBVankI expression causes the biocl( of the transition from larval to pupal stage. (A) Light nnicrographs of yw;h-Gal4 larva and 
pupa (control) and h-Gal4>TnBVankl larvae at different days AED. The scale bar is 500 ^m. (B) Larval length of different genotypes, at 96 h AED. Five 
larvae of each genotype were analyzed and as control we measured larval length of yw and h-Gal4 and UASp-TnBVankl stocks and yw;h-Gal4. Graph 
represents mean ± standard deviation (s.d.); there is no significant (NS) length difference between h-Gal4>TnBVankl (2680±83 ^m) and yw;h-Gal4 
larvae (2580±82 |im). (C) Larval length of h-Gal4>TnBVankl increases during the extended larval life. Five h-Gal4>TnBVankl larvae were measured at 
different days AED; values are the mean ± s.d. of three independent experiments. The mean values of h-Gal4>TnBVankl larval length at four and 
eighteen days AED are shown above the bars. 
doi:1 0.1 371 /journal.pone.00951 04.g001 



ethanol persisted as third instar (n = 30). This result confirms that 
the developmental arrest of TnBVankI larvae is due to a reduced 
level of 20E. However, the rescued pupae failed to develop into 
adult flies. This may be due to the fact that the large peak of 20E 
required to trigger metamorphosis is not generated by TnBVankI 
pupae and cannot obviously be supplied with food at this 
developmental stage. 

It has been reported that a positive feedback is required for the 
transcriptional up-regulation of enzymes acting at late steps in the 
ecdysone biosynthetic pathway [28]. Therefore we analyzed the 
expression of Disembodied (Dib), the downstream step enzyme 
C22 hydroxylase, which appeared strongly reduced in all 
TnBVankI PGs analyzed (n = 60) compared to control 
(Figure 2E,F). This data is in agreement with low levels of 20E 
detected in TiBVankl lai-vae. 



rnBVANKI affects PG morphology 

Using a polyclonal antibody raised against two synthetic 
peptides of TnBVANKl [16], we detected the distribution of 
TnBVANKl protein in TnBVankI PGs of five days AED larvae. As 
shown in Figure 3A-C, the protein was strongly expressed and 
present only in the cytoplasm of PG cells, confined to stroke- 
shaped particles. We next analyzed the TnBVankI PG gross 
morphology. To visualize the PG we used phm-Gal4,UAS- 
mCDSr.GFP stock. PGs from control larvae (Figure 3D) were 
significandy larger (n = 50; t = 50.41; p<0.0001) (Figure 3F) than 
TnBVankI PGs (Figure 3E). In addition, the TnBVankI PG cells 
showed a cytoplasmic rather than the expected membrane 
distribution of mCD8::GFP (Figure 3C,E) [29]. Measurements of 
the PG cell area did not show any reduction induced by TnBVankI 
expression (n = 50; t= 1.262; NS) (Figure 3G). Therefore, the 
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Figure 2. TKie expression of TnBVankI in protlioracic gland affects the E biosynthesis. (A) Ring gland includes the prothoracic gland (PG; 
yellow), the corpora allata (CA; orange) and the corpora cardiaca (CC; red). (B) The expression of the TnBVankI gene is driven in the different ring 
gland compartments, highlighted in green, by three Ga/4 drivers. P0206-Gal4>TnBVankl , expressed in PG and CA, causes the developmental arrest at 
the last larval stage; aug21 -Gal4>TnBVankl (CA) does not induce any developmental defects; phm-Gal4>TnBVankl (PG) blocks the transition from 
larval to pupal stage. (C) Total 20E titer in five larvae of UASp-TnBVankl stock (white bars), phm-Gal4/TM6B (grey bars) and phm-Gal4>TnBVankl 
(black bars), at different time (hours AED). In the control stocks UASp-TnBVankl and phm-Gal4/TM6B, the 20E peak which induces the pupariation is 
present at 120 h AED. Instead, this peak is absent in phm>TnBVankl larvae at 120 h AED and during the extended larval life. Error bars represent s.d.; 
*** = p<0.0001 versus controls (UASp-TnBVankl and phm-Gal4/TM6B). The mean values of total 20E at 120 h AED of different genotype larvae are 
shown above the bars. (D) Feeding TnBVankI larvae with medium supplemented with 20E induces the pupariation (red), while TnBVankI larvae fed 
with medium containing ethanol (EtOH) do not reach the pupal stage (green). Values are the mean ± s.d. of three independent experiments. The 
yw;phm-Gal4 larvae serve as background control (blue). Immunostaining with anti-Dib in yw;phm-Gal4 (E) and TnBVankI (F) PG reveals that the 
expression of Dib is strongly reduced in all TnBVankI PGs analyzed. Panels E,F are at the same magnification and the reference scale bar is 25 |j,m 
indicated in E. 

doi:1 0.1 371/journal.pone.00951 04.g002 



observed size difference of PG can be attributed to a reduction of 
the cell number. We then assayed if apoptosis occurs, using 
Cleaved Caspase-3 antibody [30] and TUNEL labeling [31]. The 
Caspase-3 activity (Figure 3IJ; n = 60) and TUNEL positive 
staining (Figure 3L,M; n = 60) found in a few cells of TnBVankI 
PGs suggested that the occurrence of cell death during develop- 
ment can partly account for this difference, which could be related 
to the developmental arrest induced by TnBVANKl. However, 
the possibility that this protein can also disrupt PG activity cannot 
be ruled out. Therefore, to assess the relative contribution of these 
two effects, not mutually exclusive, we expressed TnBVankI in PG 
cells at different time points during larval life, using a temperature 
sensitive form of the Gal4 repressor Gal80, GalSO''" [32], that 
allows to regulate the phm-Gal4 activity. TnBVankI and control 
larvae were initially raised at 21°C, and then shifted to the 
restrictive temperature (31°C) at specific time points (96 h, 72 h 
and 48 h AED) to promote Gal4 activity. The temperature shift 
did not affect the proper development of the control larvae, which 
pupate normally. Conversely, the larvae expressing TnBVankI 



failed to pupate, increased their size and survived for an extended 
period. 

For each time point we also analyzed the PG size at 1 20 h AED 
(Figure 3N). When the TnBVankI expression was triggered at 96 h 
or 72 h, the PGs size was not significantly different from controls 
(respectively n= 10; t = 0.07636; NS and n= 10; t= 1.336; NS). 
While the earlier induction of the transgene expression, at 48 h 
AED, strongly affected the PG size, which appeared significandy 
reduced (n= 10; t= 11.68; p<0.0001). In addition, we examined 
whether by inhibiting apoptosis with ectopic expression of p35 
[33] it would be possible to rescue the phenotype produced by the 
expression of TnBVankI in the PG. Coexpression of UAS-p35 and 
UASp-TnBVankl in the same PG cells with phm-Gal4 driver did not 
rescue the developmental arrest phenotype (n = 58). Collectively, 
these data indicate that the developmental arrest induced by 
TnBVankI does not depend on the reduced PG size triggered by 
apoptosis, but on its capacity to disrupt PG functioning when 
expressed before the production of the 20E peak. 
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Figure 3. TnBVANKI distribution in the PG cells and its effects on PG. The immunolocalization of InBVANKI in PG cells (marked with 
mCD8::GFP, green), analyzed with anti-InBVANKI (cyan), shows its presence in stroke-shaped particles (A-C), which are distributed only in cytoplasm 
(nucleus is stained with Propidium Iodide, red) (B,C). (D) At five days AED, the PG of yw;phm-Gal4 larvae, marked with GFP, is significantly larger (+ 
54%) than the TnBVankI PG (E). (F) The graph represents the mean ± s.d.; 50 PGs were analyzed; *** = p<0.0001. (G) Measurement of PG cell area 
shows no difference between yw;phm-Gal4 and TnBVankI. 50 PGs were analyzed; NS: not significant. (H-J) Immunostaining with anti-Cleaved 
Caspase-3 (red) or TUNEL labeling (red) in PG cells marked with GFP. In the control yw;phm-Gal4 no caspase or TUNEL signal is detected (H,K), while in 
TnBVankI PG few cells undergo apoptosis (l,J,L,M). PGs in panels A,D,E,H,I,K,L are at the same magnification and their scale bar is 25 jim and is 
showed in A. Scale bar in B and C is 5 |im and showed in B. Boxed regions are magnified in J and M and the reference scale bar is in B. (N) Larvae of 
yw;tub-Gal80"/+;phm-Gal4/+ {Gal80"-phm-Gal4) and UASp-TnBVankl;UASp-TnBVankl/tub-Gal80";phm-Gal4/+ {Gal80"-TnBVankl) were raised at 21 C 
(cyan) for different time intervals, then shifted at 31 C (red) and dissected at 120 h AED. PG size from larvae incubated at 2T C until 96 h AED or until 
72 h AED shows no significant (NS) differences from control larvae. PG size is strongly reduced in Gal80"-TnBVankl larvae incubated at 21 "C until 48 h 
AED compared to PG from control larvae (*** = p<0.0001). Graph represents mean ± s.d.; 10 PGs were analyzed for each experiment. 
doi:1 0.1 371 /journal.pone.00951 04.g003 



7nBVANK1 affects the cytoskeletal network in the PG 
cells 

The altered TnBVankI PG cell morphology and the associated 
mislocalization of mCD8::GFP prompted us to analyze the 
cytoskeletal network in these cells. 

We investigated F-actin and 0(- tubulin distribution in TnBVankI 
PGs and we observed an altered cytoskeletal organization in all 
analyzed glands (n = 60). As shown by phaUoidin staining 
(Figure 4A-D), cortical actin did not appear regularly distributed 
in TnBVankI PG cells, in which thick masses of actin filaments 
were detected (Figure 4C,D). The microtubule network was 
investigated by analyzing the distribution of a-tubulin-GFP fusion 
protein, which was coexpressed with TnBVankI in the PG. 
Compared to control, expressing only a-tubulin-GFP protein 
(Figure 4E,F), the cytoskeleton of the TnBVankI PG cells appeared 
strongly affected, as shown by the formation of thick bundles of 
microtubules (Figure 4G,H). The dynamic function of the 
microtubule network was then analyzed in TnBVankI PGs 



(n = 60) by assessing the distribution of the minus-end-directed 
microtubule motor dynein, using an anti-Dynein heavy chain 
antibody [20]. Compared to the control (Figure 41 J), cells of 
TnBVankI PG displayed a reduced cortical distribution of dynein, 
along with some large dynein dots (Figure 4K,L). These data 
indicate that the whole cytoskeletal network is markedly altered in 
the PG cells expressing TwBVANKl. We also analyzed GalStf'- 
TnBVankl PG ceUs at different time points (96 h AED, 72 h AED 
and 48 h AED) and we observed that F-actin organization is 
strongly altered when larvae were shifted to restrictive temperature 
at 48 h AED. This suggests that the prolonged expression of 
TnBVankI during development is causative of the disruption of 
cytoskeleton (Figure SI). Moreover, as discussed above, no adult 
phenotypic effect or developmental delay was produced by the 
expression of TnBVankI in different tissues, using a wide range of 
tissue specific Gal4 drivers (Table SI). This suggests that the 
cytoskeletal structure is not affected in all tissues, as can be 
observed in the fat body (Figure S2). This is further corroborated 
by our previous study showing that in the Drosophila oocyte 
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Figure 4. TnBVankI PG cells have an altered cytoskeleton. Phalloidin staining in control (A,B) and in TnBVankI (C,D) PG cells. F-actin shows an 
altered distribution, characterized by thick masses of filaments in TnBVankI PG cells. (E-H) cz-tubulin-GFP fusion protein was expressed in yw;phm- 
Gal4 and TnBVankI PG to investigate the microtubule network. Compared to control (E,F), in TnBVankI the microtubule cytoskeleton is strongly 
affected and forms bundles (G,H). (I-L) Immunostaining with anti-Dynein heavy chain shows that, compared to control (l,J), in TnBVankI PG cells the 
cortical localization of this protein is reduced and characterized by an evident dotted distribution (K,L). For each immunostaining we analyzed 60 PGs 
of five days AED larvae. PGs in panels A,C,E,G,I,K are at the same magnification and the reference scale bar is 25 |im and showed in A. Boxed regions 
are magnified in B,D,F,H,J,I and the reference scale bar 5 jim is indicated in B. 
doi:1 0.1 371/journal.pone.0095104.g004 



TnBVANKl interferes with proper microtubule and microtubule- 
motor protein functions [16], and does not affect the overall 
cytoskeletal structure. 

TnBVANKI expression alters the cholesterol trafficking 
endocytic pathway of PG cells 

The observed negative impact of TreBVANKl on the cytoskel- 
eton of PG cells may reduce the level of ecdysteroid biosynthesis 
by disrupting the uptake, transport and trafficking of sterols, 
essential steps for ecdysteroid biosynthesis [34] . Cholesterol, which 
cannot be synthesized by insects [35], enters the steroidogenic cells 
through a receptor-mediated low-density lipoprotein (LDL) 
endocytic pathway [36], which targets cholesterol to the 
endosomes. It is then transformed into 7-dehydrocholesterol in 
endoplasmic reticulum and transported to other subi:ellular 
compartments through further metabolic steps of the ecdyster- 
oidogenic pathway [35]. We analyzed lipid vesicular internaliza- 
tion and trafiicking in the TnBVankI PG cells with a staining 
procedure using Oil Red O. Conversely to control (Figure 5A), in 
all TnBVankI PGs analyzed (n = 60), we observed a var)'ing level of 
evident increased accumulation of lipid droplets (Figure 5B,C). 
Then, using filipin, which specifically stains non-esterified sterols 
[37], compared to control (Figure 5D), TnBVankI PGs (n = 60) 
showed a marked cholesterol accumulation in discrete vesicular 
drops (Figure 5E,F). These data suggest that TwBVANKl does not 
affect lipid uptake, but that the endocytic pathway is in some way 
disrupted. 

The endocytic pathway is organized into three major compart- 
ments, each characterized by specific Rab GTPase proteins that 
can be used as tags for the different endosomes [38]. Early 
endosomes are enriched in Rab5, late endosomes are associated 
with Rab7, and Rab 11 marks the recycling endosomes. We used 
antibodies directed against these Rab proteins to investigate the 



endocytic pathway in PG cells (n = 60 PGs for ea[:h experiment) 
[19]. The cellular distribution of the early (Figure 6A,B) and 
recycling endosomes (Figure 6C,D) appeared to be comparable 
between PGs of control (Figure 6A,C) and of TnBVankI larvae 
(Figure 6B,D). Whereas, compared to control (Figure 6E), in 
TnBVankI PGs few Rab7 positive vesicles were detected 



yw;phm-Gal4 phm-Gal4>TnBVank1 




Figure 5. TnBVankI PG cells show lipids accumulation. (A) In the 

control yw;phm-Gal4 there are few lipid droplets stained with Oil Red 0, 
while in TnBVankI cells several lipid droplets are detected (B,C). (E,F) In 
TnBVankI there is also a sterol accumulation, shown by filipin staining, 
which is absent in control PG (D). 60 PGs were stained for each 
experiment. PGs in panels A,B,D,E are at the same magnification and the 
reference scale bar 50 |im is showed in A. Boxed regions are magnified 
in C,F and the reference scale bar is 5 ^m indicated in C. 
dol:l 0.1 371/journal.pone.0095104.g005 
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(Figure 6F). This suggests that T/zBVANKl may somehow affect 
the endocytic pathway. 

We then analyzed the PG distribution of endosomes carrying 
the Hepatocyte growth factor-regulated tyrosine substrate (Hrs) 
(Figure 7A). This protein regulates inward budding of endosome 
membrane and multivesicular bodies (MVBs)/late endosome 
formation [22]. Interestingly, quite a few Hrs marked vesicles in 
the TnBVankI PG cells showed the stroke-shaped form associated 
with TnBVANKl signals (Figure 7D,E). In addition, most of the 
immunodetection signals of TizBVANKl (Figure 7E) colocalized 
with the Hrs marked vesicles (Figure 7F; Pearson's coeffi- 
cient = 0.96±0. 06). In contrast, most of these vesicles showing a 
normal round shape did not colocalize with TnBVANKl. This 
finding suggests an interaction of TjzBVANKl with endosome 
associated proteins, which may partly account for the observed 
alterations of the endocytic trafficking routes. 

MVBs formation is controlled by a set of proteins, the 
endosomal sorting complex required for transport, ESCRT-0 to 
III, which sequentially associate on the cytosolic surface of 
endosomes [39]. A partner of the ESCRT proteins, which also 
regulates the making of MVBs, is the ALG-2-interacting protein X 
(Alix), first characterized as an interactor of apoptosis-linked gene 



protein 2 (ALG-2) [40]. It has been reported that the late 
endosomal lipid lysobisphosphatidic acid (LBPA) and its partner 
protein Alix play a direct role in cholesterol export [41]. 
Therefore, by using an antibody directed against Alix, we analyzed 
the distribution of this protein in control and TnBVankI PGs 
(Figure 7B,G). According to its multifunctional activity [42], Alix 
was found widely distributed in the cytoplasm of wild type cells 
(Figure 7B), and, as expected, marked some Hrs positive vesicles 
(Figure 7C). Interestingly, in the TnBVankI PG cells the 
TNBVANKI positive stroke-shaped structures showed a strong 
colocaHzation with Alix (Figure 7H; Pearson's coeffi- 
cient = 0.99±0. 07). In addition, several of these Alrx positive 
stroke-shaped structures colocalized with Hrs (Figure 71; Pearson's 
coefficient = 0.95±0.16), indicating that these are modified 
endocytic vesicles. This strong interaction of TnBVANKl with 
Alix containing vesicles and the altered cholesterol distribution 
observed in PG are concurrent evidences indicating that the 
cholesterol route was altered. Therefore, the interaction between 
TjzBVANKl and endosomes specifically affects the endosomal 
trafficking of sterols, likely limiting their supply to subcellular 
compartments where ecdysteroid biosynthesis takes place [35]. 

Discussion 
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Figure 6. TnBVANKI disrupts the endocytic patKiway in PG 
cells. 60 PGs stained for RabS (A,B), Rabll (C,D) and Rab7 (E,F) in 
yw;phm-Gal4 and TnBVankI larvae at five days AED. The distribution of 
endosomes marked with RabS (A,B) and Rabi 1 (C,D) is not affected by 
TnBVankI expression, while a reduction in number was observed for 
late endosomes marked with Rab7 (E,F). All panels are at the same 
magnification and the reference scale bar is 5 p,m showed in A. 
doi:1 0.1 371/journal.pone.00951 04.g006 



PDVs are among the major host regulation factors used by 
parasitic wasps to subdue their hosts, which show immunosup- 
pression and a number of developmental and reproductive 
alterations associated with disruption of their endocrine balance 
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Figure 7. TnBVANKI protein colocalizes with Hrs- and Alix- 
positive vesicles. Confocal images of PG of yw;phm-Gal4 (A-C) and 
TnBVankI (D-l) larvae stained for Hrs (cyan), Alix (red) and FnBVANKI 
(green). In the control cells Alix (B) and Hrs (A) are widely distributed in 
the cytoplasm and their signals partially overlap (C). In TnBVankI cells 
(D,F) a number of vesicles marked by Hrs have different shape 
compared to those present in controls (A,C). These modified vesicles 
show a strong colocalization with TnBVANKI signal (E,F), demonstrating 
that TnBVANKI protein is associated with Hrs-marked vesicles. In 
TnBVankI cells (G) most of Alix-marked vesicles have different shape 
compared to those present in controls (B). Immunostaining with anti- 
Alix and anti-TnBVANKI shows a strong colocalization of TnBVANKI and 
Alix signals in the stroke-shaped vesicles (H). In these modified vesicles 
the Alix and the Hrs signals are both detected (I). PGs in all panels are at 
the same magnification and the reference scale bar is 5 |xm indicated in 
A. 

doi:10.1371/journal.pone.0095104.g007 
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Figure 8. Schematic overview showing the step in which ThBVANKI affects the ecdysone biosynthesis in the PG cells. EE: early 
endosomes, MVB: multivesicular bodies. 
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[3,4,6] . Relatively more studies iiave addressed the host immuno- 
suppression mechanisms, focusing on virulence factors in the ank 
gene family largely shared among difiFerent taxa [43]. While an 
immunosuppressive function has demonstrated for the PDV ank 
gene family, if and how these viral genes impact endocrine 
pathways or other targets has not yet been addressed [11]. Here 
we report experimental evidence demonstrating the role of a 
TnBV ank gene in the disruption of E biosynthesis and the 
induction of developmental arrest. 

The proteins encoded by PDV ank genes show significant 
sequence similarity with members of the IkB protein family 
involved in the control of NF-kB signaling pathways in insects and 
vertebrates [44] . Because they lack the N- and C-terminal domains 
controlling their signal-induced and basal degradation, they are 
able to bind NF-kB and prevent its entry into the nucleus to 
activate the transcription of genes under kB promoters [15,45,46]. 
The ank gene family is one of the most widely distributed in PDVs 
and contains members which are rather conserved across viral 
isolates associated with difierent wasp species [10,15,46-48]. 
These genes hkely originate from horizontal gene transfer from a 
eukaryote, which could be the wasp, the host or another organism. 
Indeed, the nudiviruses, ancestors of bracovrruses [7], do not 
encode any gene showing similarity with ank family members. 
Their multiple acquisition and stabilization in different evolution- 
ary Kneages are clearly indicative of the key role they play in 
successful parasitism. This also suggests that ank genes may be 
involved in multiple tasks on host parasitization, by influencing 
different physiological pathways. 

Here, we provide experimental data that corroborate this 
hypothesis for TnBVankI, a gene of the bracovirus associated with 
the wasp T. nigriceps (TnBV), which parasitizes the larval stages of 
the tobacco budworm, H. virescem. Using Drosophila as a model 
system, we show that the TnBVANKl protein acts as a virulence 
factor disrupting E biosynthesis (Figure 8) and causes develop- 
mental arrest of the larvae, which fail to pupariate. The number of 
late endosomes is reduced in the TnBVankI expressing cells and 
this is concurrent with an interesting change of Hrs-T^BVANKl 
positive vesicle morphology. This defective mechanism in MVB 
and late endosome formation is accompanied by an evident 



alteration of sterol trafficking as indicated by the accumulation of 
lipid and sterol-rich vesicles. 

Cholesterol is processed to free cholesterol by lipase in the 
endosomal compartment and after that it moves to other 
compartments entering the ecdysone biosynthesis machinery 
[34]. Recent evidences from mammalian cell studies indicate that 
the late endosomal lipid LBPA and its partner Alix play a role in 
controlling the cholesterol export from endosomes [41]. Our 
finding that TnBVANKl interacts with Alix positive vesicles and 
affects the sterol delivery suggests that Alix function in cholesterol 
export is conserved between Drosophila and mammals. 

Our data let us to hypothesize that in TnBVankI expressing PG 
cells cholesterol may be trapped into the MVBs. This block leads 
to insufficient sterol supply to reach the ecdysone level necessary to 
complete development. Interestingly, the fact that TnBVankI 
expression in other tissues did not alter development suggests that 
TnBVANKl impact on cholesterol trafficking may deeply affect 
the PG cells engaged in an intense steroidogenic activity. 

We show that TkBVANKI disrupts the cytoskeletal structure of 
PG cells, and this appears to be a PG specific alteration. Indeed, in 
our previous work we demonstrated that the targeted expression of 
this ank gene in Drosophila germ cells alters microtubule network 
function in the oocyte, as shown by the mislocalization of several 
maternal clues, without affecting the cytoskeletal structure [16]. 
Therefore, we cannot exclude that the specific targeted effect of 
TnBVANKl on the cytoskeleton function of PG cells may have a 
negative impact on ecdysteroidogenesis. However, it can also be 
true that the disruption of the cytoskeletal structure of these cells 
could be a downstream consequence of the impaired steroidogenic 
activity. The altered cell physiology and the consequent accumu- 
lation of lipids and sterols may have wide-ranging and more 
generalized effects on cell architecture/dynamics and survival. In 
fact, the prolonged expression of TnBVankI by phm-Gal4 during 
larval development causes cytoskeleton alteration and also 
apoptosis of a few cells, which may partly account for the 
observed reduction of the PG size. 

It is interesting to note that the developmental arrest at L3 larval 
stage induced by TnBVankI expression in the PG perfectly mimics 
the developmental alteration of parasitized tobacco budworm 
larvae, which can regularly undergo larval molting but ultimately 
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fail to pupate [13,49]. The reduced gland size observed in 
parasitized larx^ae and the low basal production of ecdysteroids 
[14,50] are fully compatible with a general reduction of the 
biosynthetic activity likely induced by ank genes. However, in 
naturally parasitized larvae these symptoms are also associated 
with a disruption of PTTH signaling, which requires active T»B V 
infection of PG, where different viral genes are expressed [13,51]. 
The high similarity of the recorded phenotypes represents a solid 
background on which to design specific experiments on the natural 
host. Indeed, the results reported here set the stage for specific in 
vivo studies in parasitized host larvae, that will have to address the 
respective roles of difiFerent TnBV genes in the suppression of 
ecdysteroidogenesis . 

Supporting Information 

Figure SI Prolonged expression of TnBVankI in PG 
cells during development alters cytoskeleton structure. 

Phalloidin staining in PGs from Gal80'^-phm-Gal4 and GalSf/'- 
TuBVankl lar\'ae raised at 2 1 °C (cyan) for different time intervals, 
then shifted at 31°C (red) and dissected at 120 h AED. PG ceU 
cytoskeleton from Gal8(f -TnBVankI larvae incubated at 21°C until 
96 h AED (D) or until 72 h AED (E) shows no significant differences 
from Gal80^-phm-Gal4 (A,B). F-actin cytoskeleton is completely 
altered in PG cells of Gal8(/'-TnBVankl larvae incubated at 21°C 
until 48 h AED (F) compared to the control treated in the same 
condition (C). PG cells in all panels are at the same magnification 
and the reference scale bar 5 |im is indicated in A. 
(TIF) 

References 

1. Quicke DL (1997) Parasitic wasps. Chapman and Hall, London. 

2. Beckagc NE, Gclman DB (2004) Wasp parasitoid disruption of host 
development: impUeations tor new biologically based strategies for insect 
control. Annu Rev Entomol 49: 299-,S30. 

3. Pennacehio E, Strand MR (2006) Evolution of de\'elopmental strategies in 
parasitic hymenoptera. Annu Rev Entomol 51: 233-258. 

4. Webb BA, Beckage NE, Hayakawa Y, Krell PJ, Lanzrein B, et al. (2000) Family 
Polydnaviridae. In Virus Taxonomy; Seventh Report of the International 
Committee on Taxonomy of Viruses: 253-260. 

5. Strand MR (2010) Polydnaviruses. In: Insect Virology, Asgari, S. and K.N. . 
Johnson (eds.). 216—241. Academic Press, Norwich, UK. 

6. Webb BA, Strand MR (2005) The biology and genomics of polydnaviruses. In 
Comprehensive Molecular Insect Science 6: 323-360. 

7. Bezier A, Annahcim M, Herbinirre J, Wetterwald C., (iyapay G, et al. (2009) 
Polydnaviruses of braconid wasps derive from an ancestral nudivirus. Science 
323: 926-930. 

8. Drczen JM, Herniou EA, Bezier A (2012) Evolutionary progenitors of 
bracoviruses. In: Beckage NE, Drezen JM (eds) Parasitoid Viruses - Symbionts 
and Pathogens, Elsevier Inc: 15—31. 

9. VolkofTAN, Drezen JM, Cusson M, Webb HA (2012) The organization of genes 
encoding ichnovirus structural protein.s. In: Beckage NE, Drezen JM (eds) 
Parasitoid Viruses - Symbionts and Pathogens, Elsevier Inc: 33—45. 

10. Strand MR (2012) Polydnavirus gene expression profiling: what we know now. 
In Beckage NE, Drezen JM (eds) Parasitoid Viruses - Symbionts and Pathogens, 
Elsevier Inc: 139-147. 

11. Gueguen G, Kalamarz ME, RamroopJ, Uribe J, (lovind S (2013) Polydnaviral 
ankyrin proteins aid parasitic wasp survival by coordinate and selective 
inhibition of hematopoietic and immune NE -kappa B signaling in insect hosts. 
PLoS Pathog 9: el003580. 

12. Pruijssers AJ, FalabeUa P, Eum JH, Pennacehio F, Brown MR, et al. (2009) 
Infection by a symbiotic polydnavirus induces wasting and inhibits metamor- 
phosis of the moth Pseudoplusia includem. ] Exp Biol 212: 2998—3006. 

13. Pennacehio F, Malva C, Vinson S (2001) Regulation of host endocrine .system by 
the endophagous braconid Cardiochilcs tiigdcep.^ and its polydnavirus. In: Edwards 
JP, Weaver RJ (eds) Endocrine interactions of insect parasites and pathogens. 
BIOS Scientific, Oxford. 123-132. 

14. Pennacehio F, FalabeUa P, Sordetti R, Varricchio P, Malva G, et al. (1998) 
Prothoracic gland inactivation in Heliothis virescens (F.) (Lepidoptera:Noetuidae) 
larvae parasitized by Cardiochihs nigriceps Viereck (Hymenoptera:Braconidae). 
J Insect Physiol 44: 845-857. 

15. FalabeUa P, Varricchio P, Provost B, Espagne E, Eerraresc R, et al. (2007) 
Characterization of the IkappaB-lUce gene famUy in polydnaviruses associated 
with wasps belonging to difierent Braconid subfamihes. J Gen Virol 88: 92—104. 



Figure S2 Expression of TnBVankI in fat bodies does 
not affect cell morphology. Phalloidin staining in fat bodies 

from die controlj-ro; Isp2-Gal4; UAS-mCD8::GFP {A,B) and from fat 
bodies expressing TnBVankI Isp2-Gal4, UAS-mCDSr.GFP/ 
TnBVankI (C,D). Fat bodies are at the same magnification in all 
panels and the scale bar is indicated in A. 

(TIF) 

Table SI Effects of TnBVankI expression using differ- 
ent Gal4 drivers. 

(PDF) 
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